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摘要 : 本 文采 用 分 子 动力 学 模拟 方法 研究 离子 的 存在 对 水 茜 气 核 化 过 程 的 影响 。 模 拟 系 综 初 始 温度 为 600 K， 达 到 平衡 后 撤去 
水 蒸气 热 浴 ， 将 载 气温 度 标 定 至 200 K， 与 水 蒸气 进行 热量 交换 ， 带 走 水 蒸气 核 化 释放 的 核 化 潜 热 。 研 究 结果 表明 ， 离 子 的 存 
在 能 够 促使 核 化 的 发 生 ， 但 离子 的 浓度 对 系 综 内 最 大 团 徐 的 生长 影响 较 小 。 根 据 团 簇 分 析 统 计 系 综 内 含 20 个 NavVCL 时 ， 核 化 
率 约 为 1.01x1023/(cm3.s)， 与 经 典 核 化 理论 计算 值 4.89x10”3/(cm3.s) 在 同一 个 数量 级 。 
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Molecular Dynamic Study of the Effect of Ions on Water Vapor Nucleation 
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Abstract: Molecular dynamics simulation was employed to explore the mechanism of heterogeneous nucleation of 
water vapor under the effect of ions. After the System reached to the equilibrium state at initial temperature 600 K, 
the temperature of carrier gas was rescale at 200 K, which worked as the heat resource for the water vapor 
nucleation. The nucleation heat was removed by carrier gases during the nucleation process. The results show that, 
ions promote the nucleation progress on the initial stage, and the concentration of ions has less effect on the growth 
of the largest cluster in the system. The nucleation rate obtained 位 om the cluster analysis was 1.01x1028/(cms3.s) in 
the system with 20 Na* and 20 CL ions, which is in the same order of magnitude with the value 4.89x1028/(cma3.S) 
predicted by classic nucleation theory analysis. 
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1 前 千 气 中 悬浮 的 离子 又 分 为 小 离子 与 大 离子 两 种 ， 小 离 

子 是 几 个 中 性 分 子 团聚 在 一 个 带电 离子 周围 形成 的 
相 变 是 自然 界 广泛 存在 的 一 种 物理 现象 ， 在 能 粒子 ， 小 粒子 吸附 在 较 大 的 中 性 气 溶 胶 粒子 上 形成 
源 利 用 、 强 化 传 热 等 领域 起 着 重要 作用 02。 在 大 气 ”大 离子 。 同 一 符号 的 大 小 离子 的 浓度 在 海洋 乡村 城 
层 中 ， 当 水 蒸气 的 相对 饱和 度 达到 过 饱和 或 大 气 中 市 不 尽 相 同 ， 海 洋 上 空 小 离子 较 多 ， 城 市 上 空 大 离 
有 足够 多 的 灰 侍 粒子 时 ， 就 会 发 生 水 蒸气 凝结 过 程 子 较 多 ， 并 且 大 小 离子 浓度 变化 趋势 相反 。 海 陆 上 


而 形成 云雾 。 随 着 经 济 的 发 展 、 科 技 水平 的 提高 ， 空 的 正 负 离子 浓度 相差 不 大 ， 但 正 离子 稍 多 1]。 
气 溶 肌 颗粒物 污染 越 来 越 引起 人 们 的 注意 ， 大 气 气 Winkler 等 由 实验 研究 有 机 蒸气 在 大 小 为 1-24 
溶胶 严重 的 影响 了 空气 质量 、 危 害 了 人 类 健康 。 大 nm “的 小 团 徐 和 纳米 颗粒 上 的 异 质 核 化 过 程 并 与 经 
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典 核 化 理论 及 Kelvin 关系 对 比 ， 结 果 发 现 ， 异 质 核 
越 小 越 能 激发 核 化 。Inci 等 回 通过 引入 与 核 化 蒸气 

分 子 大 小 相同 的 粒子 作为 异 质 核 来 研究 异 质 核 化 ， 

并 改变 异 质 核 与 蒸气 分 子 的 相互 作用 力 来 讨论 其 对 
核 化 的 影响 。 研 究 结 果 表 明 ， 异 质 核 化 与 蒸气 分 子 
相互 作用 力 较 小 的 时 候 ， 异 质 核 对 小 团 簇 胚 核 的 形 
成 功 影响 很 小 ; 随 着 相互 作用 力 的 增加 ， 团 簇 形成 
自由 能 降低 。 张 超 等 采用 团 簇 分 析 法 分 析 大 于 不 
同 阔 值 的 团 复 数 随时 间 的 变化 情况 ， 通 过 最 小 二 乘 
法 分 别 拟 合 获取 3 种 模拟 工 况 下 的 核 化 速率 ， 对 比 
不 同 带电 量 离子 对 核 化 过 程 的 影响 。 结 果 表 明正 二 
价 离子 比 正 一 价 离子 周围 产生 的 电场 力 大 ， 对 核 化 
过 程 影响 大 。 因 此 离子 的 存在 对 水 蒸气 核 化 有 着 重 
要 的 影响 ， 而 正 负 离子 及 离子 浓度 对 水 蒸气 核 化 的 
微观 机 理 还 有 待 于 进一步 的 研究 ， 本 文 建立 分 子 动 
力学 模拟 模型 ， 研 究 离 子 影响 水 蒸气 核 化 的 微观 机 
理 及 核 化 的 微观 物理 过 程 。 


2 模拟 模型 与 计算 细节 


采用 Na+、CL 离 子 代表 正 负 粒 子 研究 其 对 水 蒸 
气 核 化 的 影响 。 图 1 为 带电 离子 存在 下 的 水 蒸气 核 
化 模拟 系统 图 ，5000 个 SPC/E 模型 的 水 蒸气 分 子 均 
匀 分 布 在 模拟 盒 29x29x29 ”nm 内 ，5000 个 毛 蒸 气 
分 子 作为 载 气 〈 未 画 出 )， 系 统 平衡 后 控制 载 气 温 
度 ， 载 气 可 带 走 水 蒸气 核 化 释放 的 显 热 和 潜 热 。 为 
使 图 示 清 晰 ， 水 分 子 用 点 表示 ，Na+、CL 离 子 用 放 
大 球体 表示 。 


hy 


图 1 存在 带电 离子 的 水 蒸气 核 化 模型 


Fig. 1 Schematic plot of water nucleation system with ions 


握 分 子 间 相互 作用 采用 经 典 的 L-] (12-6) 势 孙 
数 。 水 分 子 选用 SPC/E 模型 ， 作 用 势 函数 为 ， 
0 1 Cid 

me da (1) 


于 氧 原子 质量 较 小 ， 方 程 (1) 中 第 一 项 L-] 势 


Ooo 


$0) = 4 lO 2 一 
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函数 只 用 来 计算 氧 原子 间 的 相互 作用 力 。o6o。 和 seo 
分 别 为 氧 原子 间作 用 的 尺度 参数 和 势 阱 深度 。 后 一 
项 计算 各 个 带电 荷 原子 间 的 静电 力 ，g; 与 q; 分 别 为 
原子 i 与 j 的 电荷 ，rj 为 原子 i 与 j 间 的 距离 。&0 为 
真空 介 电 常 数 。 水 分 子 与 氢 原 子 间 的 相互 作用 力 由 
修正 的 L-J 势 函数 计算 ， 水 分 子 与 Na+、Cl 离子 的 
作用 势 函 数 取 自 文献 四， 根据 Lorentz-Berthelot 混合 
法 则 计算 其 与 水 分 子 上 的 氧 原子 作用 的 能 量 参 数 与 
尺度 参数 。 表 1 给 出 势 函 数 中 各 原子 作用 参数 及 电 
位 。 时 间 步 长 3 久 ， 截 断 半径 1.5 nm， 采 用 PPPM 
方法 对 长 程 静电 力 进行 修正 。 采 用 周期 性 边界 条 

件 ，Velocity Verlet 算法 求解 分 子 动力 学 方程 。 首 
先 ， 采 用 Nose-Hoover 热 浴 ， 模 拟 300 ps 使 得 系统 
在 600 K 温度 下 达到 平衡 ， 撤 去 热 浴 将 载 气 温度 标 
定 至 200 KK。 再 模拟 2000 ps 统计 参数 分 析 核 化 过 
程 。 


表 1.， 势 函数 中 各 原子 作用 参数 及 电荷 


Table 1 Parameters and charges of atoms for potential 


function 
Atom o/A 2 /kJ moll qi/e 
O 3.166 0.6502 -0.8476 
H 0 0 0.4238 
Ar 3.405 0.996 0.0 
Na 2.583 0.4184 1.0 
Cl 4.400 0.4184 -1.0 


3 模拟 结果 及 分 析 
3. 1 核 化 过 程 热力 学 参数 

为 讨论 离子 浓度 对 水 燕 气 核 化 的 影响 ， 系 综 内 
分 别 引入 20、50、100、200 个 NaVCL 离 子 对 。 图 
2 给 出 了 不 同 离子 浓度 系 综 内 水 蒸气 的 动能 与 势能 
随时 间 的 变化 关系 。 系 综 在 模拟 达到 平衡 后 撤去 热 
浴 ， 将 载 气温 度 标定 至 200 K， 在 载 气 热 浴 下 ， 水 
蒸气 温度 降低 至 核 化 平衡 温度 。 由 图 2(a) 温 度 的 分 
布 图 可 以 看 出 ， 在 核 化 初期 ， 离 子 浓度 较 低 的 系 
统 ， 核 化 慢 ， 释 放 潜 热 少 ， 温 度 下 降 略 快 ， 离子 浓 
度 较 高 的 系统 ， 系 综 内 水 蒸气 核 化 较 剧 烈 ， 核 化 释 
放 的 潜 热 未 能 被 及 时 带 走 ， 温 度 降低 略 慢 。 随 着 核 
化 进行 ， 温 度 均 趋 于 平衡 值 ， 系 综 内 离子 浓度 对 平 
衡 温度 的 影响 较 小 ， 其 对 应 的 核 化 平衡 温度 为 
372~375 ” 开 。 核 化 过 程 是 势能 降低 的 过 程 ， 由 于 离 
子 的 相互 吸引 使 得 分 子 间 距离 减 小 ， 使 得 离子 浓度 
较 大 的 系 综 对 应 的 势能 较 低 。 
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(b) 势 能 
图 2 离子 存在 时 水 蒸气 核 化 过 程 热力 学 参数 随时 间 变 化 关 
系 


Fig. 2 Evolution of thermal properties in nucleation of water 
vapor with ions (a) temperature (b) potential energy 
3.2 核 化 过 程 分 析 
Na+、CI 与 水 分 子 间 的 短程 范 德 华 力 相对 较 

弱 ， 其 主要 的 相互 作用 力 来 之 长 程 库伦 力 。 图 3 为 
系 综 内 包含 100NaVCL 正 负离子 时 水 蒸气 核 化 的 动 
态 过 程 。 与 之 对 应 图 4 为 系统 内 不 同时 刻 点 团 簇 大 
小 的 分 布 。 系 统 模 拟 300 ps 使 其 达到 600 K 高 温 的 
平衡 状态 ， 此 时 水 蒸气 分 子 均 匀 地 分 布 在 系统 内 ， 
无 稳定 团 簇 存在 。 然 后 撤去 热 浴 ， 将 载 气温 度 标定 
至 200 玉 。 随 着 载 气 与 水 蒸气 进行 热量 交换 ， 带 走 
显 热 ， 水 蒸气 系统 温度 降低 。 核 化 为 亚 稳 态 过 程 ， 
离子 与 水 分 子 的 较 强 作用 力 相当 于 一 些 扰动 ， 促 使 
核 化 发 生 ， 温 度 降低 到 核 化 点 时 ， 在 离子 周围 的 水 
分 子 以 离子 为 核 化 中 心 核 化 。 由 于 原子 间 的 相互 作 
用 随 着 距离 的 增加 而 减弱 ， 所 以 带电 离子 对 其 周围 
水 分 子 核 化 的 影响 也 随 距离 变化 。 如 图 3 中 {=500 
ps 时 水 藻 气 分 子 核 化 形成 团 簇 的 情况 ， 对 应 团 禾 大 
小 分 布 见 图 4 中 500 ps， 这 时 系统 内 大 量 小 团 簇 
形成 ， 不 稳定 的 小 团 簇 在 分 子 相互 作用 下 破裂 ， 稳 
定 的 小 团 簇 继续 长 大 。 当 t=1300 ps， 形 成 的 小 团 簇 
开始 在 相互 作用 下 案 并 形成 大 团 禾 ， 最 大 团 簇 约 由 
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500 个 水 分 子 组 成 ， 此 时 仍 有 一 些小 的 不 稳定 小 团 
徐 存 在 ， 不 断 地 破裂 又 形成 ， 系 统 内 水 燕 气 饱和 度 
降低 。 当 人 二 2000 ”ps 时 ， 系 统 内 的 核 化 基本 达到 一 
个 动态 平衡 ， 几 个 较 大 的 团 簇 稳定 存在 ， 最 大 团 簇 
约 包含 700 个 水 分 子 ， 水 分 子 单 体 均匀 分 布 在 系统 
内 。 继 续 模拟 到 2300 ps， 系 统 内 的 团 复 大 小 分 布 变 
化 较 小 ， 说 明 核 化 系统 达到 了 核 化 与 蒸发 的 动态 平 
衡 状态 。 


Em X {t=2000ps od X 


图 3 系 综 内 水 蒸气 核 化 的 动态 过 程 〈100NarWCL 离 子 ) 


Fig. 3 The snapshots for dynamic process of water vapor 


t=2300 ps 


nucleation with 100 Nar ions and 100Cl ions 
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图 4 系 综 内 团 徐 大 小 随时 间 的 分 布 〈100NaVCL 离 子 ) 


Fig. 4 Evolution of cluster size distribution in water nucleation 
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process with 100 Nar ions and 100CL ions 
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Fig. $5 Evolution of the largest cluster size distribution in water 


pa 


nucleation process with ions 

同时 讨论 了 离子 浓度 对 水 蒸气 核 化 的 影响 ， 分 
析 不 同 离子 浓度 下 ， 系 统 内 最 大 团 艇 包含 的 分 子 数 
随 核 化 时 间 的 变化 ， 如 图 5 所 示 。 在 核 化 初期 ， 离 
子 的 存在 能 够 促使 核 化 的 发 生 ， 团 复生 长 较 快 ， 图 
中 发 生 阶 路 为 小 团 簇 聚 并 到 大 团 禾 ， 最 大 团 簇 急 剧 
长 大 。 随 着 核 化 进行 ， 小 团 艇 聚 并 长 大 ， 系 统 达到 
新 的 平衡 状态 ， 最 大 团 簇 包含 约 700 水 分 子 。 不 同 
离子 浓度 的 系统 ， 核 化 达到 平衡 时 ， 最 大 团 簇 大 小 
基本 相同 ， 离 子 浓度 影响 较 小 。 
3. 3 核 化 率 计算 

核 化 率 定义 为 单位 时 间 单 位 体积 内 生产 的 活化 
核 ( 即 能 够 稳定 存在 并 长 大 的 分 子 团 簇 )。 本 文 以 包 
含 20Na*/Cl 离 子 的 核 化 系 综 为 例 ， 分 别 采用 经 典 核 
化 理论 预测 及 团 复 分布 统计 分 析 离 子 存在 时 水 蒸气 
的 核 化 率 大 小 。 核 化 系统 内 对 应 的 饱和 度 约 为 
2.9。 离 子 存 在 时 水 落 气 核 化 达到 动态 平衡 时 对 应 的 
温度 约 为 375 K， 此 温度 下 对 应 的 饱和 气 与 饱和 液 
的 密度 分 别 为 6.27x104 gcm3 和 0.9165 g:-cm3， 参 
考 表面 张力 36.94 Nm-!1。 由 于 本 研究 中 的 带电 离子 
没有 考虑 大 小 ， 采 用 均 质 核 化 的 核 化 率 计算 公式 ， 


ee AG 
= PP exp( | (2) 
B 


计算 得 出 的 核 化 率 的 大 小 为 4.89x1028/(cm3.s)。 
此 外 根据 核 化 率 定 义 ， 可 从 团 簇 分 析 方 法 求 得 核 化 
率 ， 图 6 给 出 了 包含 分 子 数 大 于 一 定 阔 值 的 团 簇 团 
簇 数 随时 间 的 变化 ， 阔 值 分 别 为 10、20、30、40、 
50。 当 阔 值 为 10 时 ， 系 综 内 团 艇 数 先 迅 速 增加 后 
降低 ， 变 化 不 稳定 。 当 阔 值 大 于 20 时 ， 团 徐 数 增 
加 的 斜率 基本 相同 ， 从 而 可 以 根据 其 变化 的 斜率 统 
计 获 得 核 化 率 。 图 中 的 虚线 为 团 复生 长 阶段 的 线性 
拟 合 ， 此 时 对 应 的 系 综 内 包含 20NaVCL 离 子 时 水 蒸 
气 核 化 的 核 化 率 约 为 1.01x1028(cm3a s)， 与 经 典 核 化 
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理论 预测 值 在 同一 数量 级 ， 高 于 文献 I 计算 得 出 的 
均 质 核 化 的 核 化 率 值 2.8x1027 /(cm3.s)， 说 明 离 子 的 
存在 对 水 蒸气 核 化 有 促进 作用 。 采 用 同样 的 方法 根 
据 团 簇 分 析 得 出 系 综 内 包含 50、100、200 ”Nat/CL 
离子 时 对 应 的 核 化 率 依次 为 ，1.54x1028/(cm3 $s)、 
2.27x1028/(cm3 s)、4.52x1028/(cm3 s)， 核 化 率 随 离子 
浓度 的 增加 而 非 线 性 的 增加 。 
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图 6 系 综 内 团 复数 随时 间 的 变化 (20NaVCL) 


Fig. 6 Evolution of clusters larger than a threshold size in the 


nucleation process with 20 Na* and 20CL ions 
4 结论 


本 文 建立 了 NaVCL 离 子 存 在 时 水 蒸气 核 化 的 分 
子 动力 学 模型 ， 模 拟 研究 NaVCL 离 子 及 离子 浓度 对 
水 蒸气 核 化 的 影响 ， 分 析 核 化 过 程 温度 、 势 能 的 变 
化 ， 根 据 团 簇 分 析 及 经 典 核 化 理论 预测 离子 存在 时 
水 蒸气 的 核 化 率 。 结 果 发 现 ， 在 核 化 初期 ， 离 子 的 
存在 能 够 促使 核 化 发 生 ， 但 随 着 核 化 进行 至 动态 平 
衡 时 ， 离 子 浓度 对 系 综 内 最 大 团 徐 大 小 影响 较 小 。 
以 含 20NarVCL 离 子 系统 为 例 ， 根 据 核 化 率 定义 由 团 
徐 分 析 统 计 得 出 核 化 率 约 为 1.01x1023/(cm3 s)， 高 于 
均 质 核 化 系统 ， 且 与 经 典 核 化 理论 计算 值 
4.89x1028/(cm3 s) 在 同一 个 数量 级 。 
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